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Abstract 

Recently we proposed a formalism, based in nuclear statistical spectroscopy, 
for efficient computation of nuclear level density, or densities of states, through 
a sum of partitioned binomial functions (SUPARB). In this Letter we extend 
the formalism to the calculation of locally averaged expectation values, with 
specific application to spin-cutoff factors and the angular momentum depen- 
dence of the nuclear density of states. 
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In a recent Letter |J we proposed a computationally efficient method to compute nuclear 
level densities, based in the ideas of nuclear statistical spectroscopy |@-f§]. While we refer 
interested readers to our original Letter for details, we briefly summarize here. Consider the 
density of states of a Hamiltonian H, 

p(E) =trd(E - H). (1) 

It is useful to partition the finite model space into subspaces, labeled by Greek letters 
a,/3, etc., each with an associated projection operator P a . Then one can define partial or 
configuration densities: 

Pa (E)=trP a 5(E-H). (2) 

The total density is the sum of the partial densities. (NB:We always include 2 J + 1 degen- 
eracies and so are formally considering state densities.) 

If the subspaces are single-particle configurations, e.g., (Og^) 4 , (0g?5/2) 2 (1si/2) 2 , etc, then 
the partial, or configuration, moments up to fourth order for any system with Z protons and 
N neutrons in the valence space can be computed directly from the one+two-body matrix 
elements of H [|5]||. For any partition (configuration) a let d a = tr P a be the dimension 
of the subspace a, and define the configuration average over the subspace to be (. . .) a = 
d a ~ trP a (. . .). Then E(a) = (H) a is the (configuration) centroid, 7(a) = ((H — E a ) 2 )]/ 2 
the configuration width, m 3 (a) = ((H — -E Q ) 3 ) a /7 3 (a) the scaled (dimensionless) third 
configuration moment, and m 4 (a) = ((H — E a ) 4 ) a /^(a) the scaled fourth configuration 
moment. 

With these moments in hand, we model the partial densities as binomial distributions, 
following a recent suggestion of Zuker 0. Starting with the binomial expansion of (1 + X) N 
and representing (^j with gamma functions, one can derives a continuous distribution, 
n / F \ _ \E x /e Y(E max /e + 1) 

P[ x) T(E x /e + l)T((E max -E x )/e + l) [6) 



where E x is the excitation energy and E max = 7(1 + X)yN/X is the maximum excitation 
energy in the binomial distribution. The binomial is appealing because one can easily 



compute the scaled third and fourth central moments: 777,3 = and m 4 = 3 — + j^, 
allowing one to control the shape of the binomial through N and A. 

To compute the density of states, we take the following steps: (1) We compute the 
configuration moments up to 3rd or 4th order. (2) We model the partial density for each 
configuration as a binomial. The binomial parameters iV and A, as well as the overall energy 
scale and centroid, are fitted to the configuration moments. (3) The partial densities are 
summed to yield the total density of states. Because of these ingredients, we will refer to our 
approach as SUPARB (SUm of PARtitioned Binomials) state densities. We have shown that 
our method models well the density of states, when compared against exact calculations, 
and that one needs third configuration moments and occasionally, but not always, fourth 
configuration moments to achieve accurate results [l] . Incidentally, none of this is special to 
the atomic nucleus. One could easily apply it to atomic electrons as well. 

We now turn to the computation of expectation values of some general operators, O, 
and define the energy- dependent locally averaged expectation value (LEV) as 

We also have configuration LEVs: 

(dm. = tcPad6{E -? ) . (5) 

X V " tr P a 5(E-H) 

Our strategy for computing the total (O(E)) is the same as for the total density of 
states. In addition to the configuration moments, we compute the weighted averages (0) a , 
(0(H — E a )) a , and {0(H — E a ) 2 ) a . These can be computed directly from one+two-body 
matrix elements similar to those for the moments P,fJ. 

The only tricky point is if O and H do not commute; where does one insert the projection 
operator P a in eqn. (|D? Because J2a Pa = 1 by the completeness of the projection operators, 
we conclude consistency is the only requirement, and take tr P a OH k , k = 0, 1, 2. 

Assume within any subspace a quadratic energy dependence, that is, 

(6(E)) = O + 1 (E - E)h + 2 (E - Efh 2 (6) 



(dropping for the moment the subspace label a). Then 



(6) = O + O 2 , (7) 
(6{H-E))= 1 (Oi + 2 m 3 ), (8) 
(6(H - Ef) = 7 2 (O + O im3 + 2 m 4 ). (9) 



Solving for coefficients, 



Oo = (Q>(m 4 ~ ml) + m 3 (Q(H - E))h - (Q(H - E) 2 )/l\ 

777,4 — ml - f 



(m 4 - l)(0(H - E))h + m 3 {(0(H - Ef)ll 2 ~ W 
u i — 2 1 ' 

Q2 = (Q(H - Ef)h 2 - m 3 (6(H - E))h ~ (g) (12) 

7774 — 7TI3 — f 

If we limit ourselves to only a linear dependence (assume O2 = 0), then 

Oo = <6>, (13) 
1 = (6(H-E))/ r (14) 

Note: in terms of the binomial parameters N, A, one finds that 7774 — m\ — 1 = 2(1 — iV -1 ). 
The generalization to the full SUPARB case, summing over partitions, is easy: 

(6(E)) = (15) 

In Figures 1-3 we illustrate our method, comparing to "exact" shell model calculations 
in full Ohu spaces. For s<i-shell nuclides we used the Wildenthal USD interaction ||, and 
compared against direct diagonalization. Fig. 1 shows the LEV of Q ■ Q and S 2 for 20 Ne, 
while Fig. 2 shows the LEV of J 2 (which will be important for the spin-cutoff factor below) 
for 22 Na, 23 Mg, and 32 S. 

Fig. 3 compares the LEV of J 2 for 48 Cr and 54 Fe in a full Ohu p/-shell calculation; here 
the "exact" calculation was through Monte Carlo sampling of path integrals |9||ll| . To avoid 



the well-known sign problem |TTJ we fitted a schematic multipole-multipole interaction to 



the T = 1 matrix elements of the FPD6 interaction of Richter et al. [pT2| . Clearly the full 



SUPARB calculation, with a quadratic dependence in each partition, works very well, and 
the quadratic is a significant improvement over the linear approximation. One could go to 
a cubic dependence, but computation of the necessary moments, while possible, would be 
very time consuming and there is not much room for improvement. 

Now we apply our formalism to spin-cutoff factors for state densities. Define the J- 
dependent density as 

p(E, J) = tr (8(E - H)6(J(J + 1) - J 2 )) . (16) 

Traditionally this is factorized to 

p(E,J)=p(E)(2J + l)Q(J,E), (17) 

and one assumes a weighted Gaussian form for Q: 

QUE) = (2J+1) ^^exp f- (J tin (18) 

We have normalized / dJQ(J, E) = 1. Here cr(E) is the "spin-cutoff" factor and is dependent 
on the energy. The spin-cutoff factor has been considered before in the context of nuclear 
statistical spectroscopy |Tj|, but the moments were computed by a random sampling of 



representative vectors and the level densities were approximated by Hermite polynomials 
(which do not guarantee nonnegative densities). 
Given the form ([TJD one finds that 

(J 2 ) = 3a 2 - \ (19) 

(The factor of 3 is because we must include the 2J+ 1 degeneracy in our traces.) Therefore, 
the problem of describing p(E, J) with SUPARB reduces to computing (J 2 ) as a function 
of energy, which were given in figures 2 and 3. 

In figure 4 we plot the exact and SUPARB J-projected level densities (without the 2 J + 1 
degeneracy) of 32 S, for J = 0, 1, 2 and 8. The results are very good. Incidentally, one could 
compute the spin-cutoff factors and thus the J-projected level densities using Monte Carlo 
evaluation of path integrals just as easily, although this has not yet been done. 



In summary, we have extended our previous SUPARB technique for state densities to 
expectation values of operators, including application to spin-cutoff factors. Our results 
look very good. Other possible application include estimating the contamination of spurious 
states in cross-shell calculations, and computation of total strengths and energy-weighted 
sum rules for transitions. 

This work was performed under the auspices of the Louisiana Board of Regents, contract 
number LEQSF(1999-02)-RD-A-06; and under the auspices of the U.S. Department of En- 
ergy through the University of California, Lawrence Livermore National Laboratory, under 
contract No. W-7405-Eng-48. 
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FIGURES 
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FIG. 1. Comparison of exact shell model locally-averaged expectation value of (a) Q ■ Q and 
(b) S 2 in 20 Ne. The histogram is from direct diagonalization, and we compare against both linear 
(dashed line) and quadratic (solid line) SUPARB estimates. 
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FIG. 2. Comparison of exact shell model locally-averaged expectation value of J 2 in (a) 22 Na, 
(b) 23 Mg, and (c) 32 S. The circles are from direct diagonalization, and we compare against both 
linear (dashed line) and quadratic (solid line) SUPARB estimates. 
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FIG. 3. Comparison of exact shell model locally-averaged expectation value of J 2 in (a) 48 Cr, 
and (b) 54 Fe. The histo gram is from Monte Carlo sampling of a path integral, and we compared 
against both linear (dashed line) and quadratic SUPARB (solid line) estimates. 
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FIG. 4. Comparison of J-projected state densities in 32 S for J = 0, 1, 2, and 8. The circles 
are 'exact' from direct diagonalization, and the solid lines are SUPARB estimates with spin-cutoff 
factor. 
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